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ABSTRACT 

This Letter explores influences of intracluster magnetic fields ( > 1/xG) submerged in the hot electron 
gas on classic Sunyaev-Zel'dovich effect (SZE) and thermal bremsstrahlung in X-ray emissions. As the 
Larmor frequency is much higher than all collision frequencies, the presence of magnetic field may lead 
to an anisotropic velocity distribution of hot electrons. For the two-temperature relativistic Maxwell- 
Boltzmann distribution, we compute modifications to the classical thermal SZE. Intracluster magnetic 
fields tend to enhance the SZE with steeper radial variations, which bear important consequences for 
cluster-based estimates of cosmological parameters. By applying the magnetic SZE theory to spectral 
observations of SZ and Chandra X-ray emissions from the galaxy cluster Abell 2163, a ~ 30 — 40/iG 
central core magnetic field Bq is predicted. For the SZ and Chandra X-ray spectral observations of the 
Coma cluster, our theoretical analysis is also consistent with an observationally inferred Bq < 10/iG. 
As the magnetic SZE is redshift z independent, this mechanism might offer a potentially important 
and unique way of probing intracluster magnetic fields in the expanding universe. 
Subject headings: cosmic microwave background — cosmology: theory — galaxies: clusters: general 
— magnetic fields — plasmas — radiation mechanisms: general 



1. INTRODUCTION 

The Sunyaev-Zel'dovich effect (SZE) in galaxy clusters 
offers a unique and powerful observational tool for cos- 
mological studies. There has been persistent progress 
in detecting and imaging the SZE in clusters. In view 
of this rapid development in SZE observations, several 
important physical effects associated with the classical 
thermal and kinetic SZEs (Sunyaev & Zel'dovich 1969, 
1980) have been further explored for their diagnostic po- 
tentials, such as relativistic effects (Rephaeli 1995), the 
shape and finite extension of a galaxy cluster with a poly- 
tropic temperature (Puy et al. 2000), halo rotation SZE 
(Cooray & Chen 2002; Chluba & Mannheim 2002), Bril- 
louin scattering (Sandoval- Villalbazo & Maartens 2001), 
early galactic winds (Majumdar et al. 2001) and cooling 
flows (Schlickeiser 1991; Majumdar et al. 2001; Koch et 
al. 2002). 

Intracluster magnetic fields have been measured using 
a variety of techniques and diagnostics, including syn- 
chrotron relics and halo radio sources within clusters, 
inverse Compton X-ray emissions from clusters, Faraday 
rotation measures of polarized radio sources either within 
or behind clusters, and cluster cold fronts in X-ray im- 
ages (Clarke et al. 2001; see Carrili & Taylor 2002 for 
a recent review). These observations reveal that most 
cluster atmospheres are substantially magnetized with 
typical field strengths of > 1/iG and with high areal fill- 
ing factors out to Mpc radii. In the cores of 'cooling flow' 
clusters (Eilek & Owen 2002; Taylor et al. 2001) and at 
cold fronts of merging clusters (Vikhlinin et al. 2001), 
magnetic fields may gain intensities of ~ 10 — 40/^G and 
thus become dynamically important. 
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Magnetic fields in the intracluster gas allows for parti- 
cle acceleration processes which modify specifics of heat- 
ing processes, such that the electron energy distribu- 
tion differs from the Maxwcll-Boltzmann distribution. 
Such stochastic acceleration processes include shocks 
and magnetohydrodynamic (MHD) waves, etc. The 
bremsstrahlung from a modified Maxwell-Boltzmann 
electron gas might account for the observed X-ray spec- 
tra up to highest energies of current X-ray observa- 
tions (Ensslin et al. 1999; Blasi 2000a). If energy in- 
jections by MHD waves are turned off, a galaxy clus- 
ter gradually thermalizes with electrons approaching a 
Maxwell-Boltzmann distribution on a rough timescale of 
~ 10 7 — 10 8 yrs. As all collision frequencies (Nicholson 
1983) are much lower than the electron Larmor frequency 
for the magnetized intracluster gas (Sarazin 1988), the 
electron velocity distribution is likely to be anisotropic as 
long as the parallel (relative to magnetic field B) pressure 
is not too much higher than the perpendicular pressure 
(Parker 1958; Hasegawa 1975). 

We presume the result of a partial electron ther- 
malization is a two-temperature relativistic Maxwell- 
Boltzmann distribution, i.e. an anisotropic velocity dis- 
tribution. This two-temperature does not mean an elec- 
tron gas having two components with different temper- 
atures, but refers to the same population with different 
average kinetic energies along and perpendicular to the 
local magnetic field. The main thrust of this Letter is 
to advance magnetic SZE theory in contexts of Chan- 
dra X-ray and radio SZE spectral observations for galaxy 
clusters. 

In Section 2, we calculate the X-ray emission and SZE 
spectra using the two-temperature relativistic Maxwell- 
Boltzmann distribution for electron velocity. Based on 
both Chandra X-ray and SZ spectral observations, we 
offer a specific prediction for the galaxy cluster A2163. 
Finally, we discuss cosmological implications of our mag- 
netic SZE theory in Section 3. 
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2. MAGNETIC SUNYAEV-ZEL'DOVICH EFFECT 

Both X-ray emission and SZE spectra are sensitive to 
the hot electron energy distribution. By the presence of 
B > 1/iG, electrons thcrmalize their parallel and perpen- 
dicular (relative to B) kinetic energies separately with a 
resulting two-temperature Maxwell-Boltzmann distribu- 
tion on timescales of ~ 10 7 — 10 8 yrs. We adopt a two- 
temperature thermal relativistic Maxwellian electron ve- 
locity distribution p c (0i,02), namely 



d/3id/3 2 



where 9i = k B T\\ / (m c c 2 ) , 82 = k B T±/(m e c 2 ), k B = 
1.38 x 10 -16 erg K _1 is the Boltzmann constant, c is the 
speed of light and m e is the electron mass; T11 and T± arc 
parallel and perpendicular temperatures, respectively; 
01 = v\\/c, fa = v±/c, 2 = 01+01 7 - (1-0 2 )- 1 ' 2 and 
7i = (1 — /3|) -1 / 2 for i = l,2 with v\\ and v± being paral- 
lel and perpendicular velocities, respectively. We assume 
B to be random over the entire cluster on scales larger 
than a typical coherence length of ~ 1 — lOkpc. Thus 
microscopically anisotropic electrons are macroscopically 
isotropic, analogous to a demagnetized ferromagnet. In- 
tegrating in all directions, the electron speed or energy 
distribution becomes 
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where TV is a normalization factor computed numerically. 

The firehose stability (Parker 1958; Hasegawa 1975) 
for velocity anisotropy requires 
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B 2 /(2ir) > n e {(mvp - (mvi)) > , 



(3) 



where n c is the electron number density, m = -fm c , 
angled brackets indicate averages over p c {0) and the 
mean temperature T is defined by 3k B T/2 = (2 (mvi) + 
(mv 2 ))/2. It then follows that 



k B T\\ : 
fc R T 



(mv\) = k B T + B 2 /(3irn c ) , 
ee (mvi) = k sT - B 2 /(6nn c ) 



with B 2 /(2irn e k B ) being the upper bound for the tem- 
perature difference Tj| — Tj_. For an observed X-ray en- 
ergy spectrum and an empirically inferred B distribu- 
tion in a cluster, one may estimate Tj| and T± by fitting 
the spectral data. By correlations between X-ray surface 
brightness and Faraday rotation measure (Dolag et al. 
2001), a power law B oc [n c (r)] a was inferred with an 
exponent a estimated from the slope of In B versus In n e 
relation. For example, one finds a ~ 0.9 for the galaxy 
cluster A119 and a ~ 0.5 for the galaxy cluster 3C 129 
with a larger uncertainty in the latter. For a = 0.5, 
T|| and T± can remain constant in a magnetized galaxy 
cluster. 

The X-ray emission rate per unit volume per unit en- 
ergy interval is given by 



jx(E x ) =n 2 c c J d0p e (0)0a B (0,Ex), 



(4) 



where a B (0,Ex) is the differential cross section (Haug 
1997) for the bremsstrahlung of an X-ray photon with 



energy Ex from an electron of speed c0. Given a cen- 
tral Bo, we fit an observed X-ray energy spectrum with 
the electron distribution (2) and Tj| — T± constrained 
by marginal firehose stabilities (Parker 1958; Hasagawa 
1975). 

To scatter a cosmic microwave background (CMB) 
photon (Birkinshaw 1999) of frequency Vi off an isotropic 
distribution of thermal electrons with speed c0, the prob- 
ability for the scattered photon of frequency i^(l + s) is 



P(s,0) = 



* (i + /V) 



x[l + Ai y 2 + (i-M 2 )(i-A/ 2 )/2]dM, (5) 

with n' = [(1 + s)(l - 0n) - l]/0, where m = (s - 
/?)/[(l + s)0] and fi2 = 1 for s > and m — —1 and 
M2 = (s + /?)/[(! + s)0] for s < 0. As the intracluster 
electron gas is macroscopically isotropic and has a thin 
optical depth r ~ 10~ 2 , integral (5) is applicable to 
the magnetic SZE analysis. For photons scattered by 
an electron distribution of expression (2), the resulting 
distribution in the fractional frequency shift s is 
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For CMB photons scattered by a hot intracluster electron 
gas, the change in the CMB spectrum at frequency v 
caused by the magnetic SZE is 
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where x = hu / (k B TcMB) , h = 6.63 x 10 _27 erg s is the 
Planck constant, Tqmb is the present CMB tempera- 
ture and r e = OT-ZVe with ctt = 6.65 x 10 -25 cm 2 being 
the electron Thomson cross section (Rybicki & Lightman 
1979) and Ac being the column density of free electrons 
along the line of sight. 

In the extensively used beta-model for intracluster hot 
electron gas (Sarazin 1988; Fabian 1994), the empirical 
radial distribution of n c (r) is 



n c (r) = n e o[l + {r/r c Y 



-3/3 c /2 
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where r c is the core radius and n e o is the central electron 
number density. It follows that t c at radius r is given by 



r e (r) =r c0 [l + (r/r c ) 2 ] 



21 (l-3/3 c )/2 



(9) 



where t o0 = n c0 aTr - cv / 7rr(3/3 c — l/2)/r(3/3 c ) and T is the 
gamma function. 4 We use eqns. (2), (5), (6) and (9) in 
integral (7) to compute magnetic SZE spectrum at r. 

For spectral observations of A2163, we estimate the 
lower limit of intracluster B from data using our mag- 
netic SZE theory. Being one of the hottest clusters, 
the thermal electron gas trapped in A2163 (rcdshift 
z = 0.203) has a mean temperature k B T c ~ 12.4±0.5keV 
and a central density n c o — 6.82 x 10~ 3 cm~ 3 with a 
core radius r c = 0.269 ± 0.025Mpc (Hubble constant 
ftioo = 0.71) and a C = 0.616 ± 0.031 (Elbaz et al. 
1995; Markevitch et al. 1996; Markevitch & Vikhlinin 
2001). By these estimates, the optical depth towards the 

4 Here, /3 C instead of j3 is used to avoid notation confusions with 
the dimensionless thermal electron speed (3. 
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Fig. 1. — The SZE spectrum of A2163 (triangles with error 
bars are data points). The heavy solid line (T = 12.4 keV) is the 
thermal SZE with parameters determined from X-ray observations 
with Bo = 0. The dash, dotted and dash-dotted lines represent 
magnetic SZE with different central magnetic field strengths Bo 
and with other cluster parameters being the same. We take a = 
0.5, so that Tj_ and Tii remain constant in the cluster. The grey 
line (T = 14.4keV, Bo = 0) coincides with the dash-dotted line. 
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Fig. 2. — Energy spectrum of the central r < 5' region of A2163 
by Chandra. The grey and dashed lines are the T = 12.4 and 
14.4 keV thermal bremsstrahlung fits. The dotted line shows our 
magnetic model of essentially the same fit as T = 12.4 keV yet with 
different Bo, Tu and T±_ of Fig. 1. The lower energy part ( < 2.1 
keV) is contaminated by the soft excess (Markevitch & Vikhlinin 
2001). The spectral resolution for energy > 5.0 keV is not high 
enough and is blended with the Fe K a line (peaked at 5.3 keV for 
z = 0.203). We fit the 2.1 - 5.0 keV band with uncertainity in T 
of 0.9 keV (90%CL). 
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Fig. 3. — Radial features of thermal and magnetic SZEs modeled 
for A2163 at u = 30GHz. The solid line is the thermal SZE, 
the dash, dotted and dash-dotted lines are models with the same 
central magnetic field strength Bo = 36^tG but different a values. 



cluster center is t c q = 0.0133 and the central Compton 
parameter is ytho = 3.21 x 10~ 4 . 

A conspicuous SZE spectrum has been observed in 
A2163 at four frequencies by BIMA at 30GHz (LaRoque 
et al. 2002), by DIABOLO at 140GHz (Desert et al. 
1998) and by SuZIE at 140, 218 and 270GHz (Holzapfel 
et al. 1997) with dust-corrections (LaRoque et al. 2002). 
These data were previously fit with a Compton param- 
eter y t h = (3.56tg;4jtg;?9) x 10" 4 for a thermal SZE 
together with a kinetic SZE for a positive peculiar veloc- 



ity V p 



= 415 



+ 1030+460 
-850 -440 



km s- 1 (68%CL) (Carlstromet al. 

2002) . A positive- velocity kinetic SZE leads to an overall 
downward shift of the thermal SZE spectrum especially 
around the zero point (~ 218 GHz) of the SZE. 

We fit the observed X-ray spectrum of Fig. 2 for 
a mean temperature T = 12.4keV by eqns. (2) and 
(4) with different B strengths and infer relations among 
B, T|| and T±. We here take a = 0.5 with Tj_ and Tp 
being constants (see parameters of Fig. 1). Inserting 
these relations into eqns. (2), (6) and (7), we obtain the 
spectra of magnetic SZE. 

Shown in Fig. 1 are SZE spectra computed for A2163 
with different parameters in comparison with the ob- 
served SZE spectrum. The model without magnetic field 
(solid line for T — 12.4keV) underestimates signals of 
SZE (especially at v = 140GHz) with x 2 = 3.82 in a 
chi-square fit. Intracluster magnetic field tends to en- 
hance signals of SZE and the best fit is Bq — 36/iG with 
X 2 = 0.78. As the marginal firehose stability (Parker 
1958) is used here, the fitting estimate represents the 
lower limit of Bq = 36/iG for A2163. Magnetic field also 
increases the null frequency of the pure thermal SZE, 
similar to the kinetic SZE with a positive velocity This 
degeneracy can be removed by SZE signals at other fre- 
quencies (e.g. 100 and 400GHz, etc.). The SZE spectrum 
data can be fit with T = 14.4keV (grey line in Fig. 1) 
and Bq — 0, but the thermal spectrum data does not fit 
well with T = 14.4keV (dashed line in Fig. 2). Shown 
in Fig. 3 are computed radial SZE features of thermal 
(B = 0) and magnetic (a ^ 0) SZEs modeled for clus- 
ter A2163 at 30GHz. With a smaller exponent a, SZE 
signals steepen from the central to peripheral parts of a 
cluster. This can be utilized to determine macroscopic 
mean B structures by obtaining spatially resolved mag- 
netic SZE intensity maps (Carlstrom et al. 2002). Note 
that 30GHz of Fig. 3 is just an illustrating example. 

3. DISCUSSIONS 

Contrary to recent results (Koch et al. 2003; Zhang 

2003) , we find that the anisotropic velocity distribution 
of electrons caused by magnetic field B enhances the 
SZE. Our model results of Figs. 1 — 3 can be critically 
tested against more precise spectral SZE measurements 
of A2163 in the frequency bands of ~ 50 — 130GHz and 
~ 300 - 600GHz by MAX, MSAM and SuZIE types of 
experiments in the frequency passbands of 90 — 670GHz 
and by AMiBA in the band 84 - 104GHz, Nobeyama 
at 21 and 43GHz, JCMT at 350 and 650GHz, SZA in 
the bands of 26 - 36GHz and 85 - 115GHz, BIMA and 
OVRO in the band of 26-36GHz, MINT at 150GHz and 
ACT at 150, 220, and 270 GHz. Multi-frequency projects 
such as the upgraded MITO (Lamagna et al. 2002; De 
Petris et al. 2002) and the OLIMPO (Masi et al. 2003) 
experiments are very promising to provide some results. 
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A2163 may involve merger shocks that could amplify 
B (e.g. Markevitch & Vikhlinin 2001). The spectral in- 
dex maps of A2163 show a spectral steepening from the 
central to peripheral radio halo regions, implying a ra- 
dial decrease of B in reacceleration models (e.g. Feretti 
et al. 2003). It was attempted to fit the SZE spectrum of 
A2163 with a combination of thermal and non-thermal 
electrons (e.g. Colafrancesco et al. 2003), but no evi- 
dence was found for hard X-ray excess due to the non- 
thermal component in the BeppoSAX data (e.g. Feretti 
et al. 2001). 

Based on X-ray and SZE measurements, 41 galaxy 
clusters were used to independently estimate Hubble 
constant ft-ioo = 0.61 ± 0.03 ± 0.18, where the uncer- 
tainties are statistical and systematic at 68% confidence 
level for Qm — 0.3 and f^A = 0.7 cosmology (Carl- 
strom et al. 2002; Reese 2003). Our analysis of A2163 
shows that intracluster magnetic field induces micro- 
scopic anisotropies in electron velocity distribution to 
enhance the SZE. It appears that inferences from clus- 
ter models without magnetic field would systematically 
underestimate ho as in the case of A2163 for which 
Holzapfel et al. (1997) inferred a lower h wo = 0.60±0.04 
against the current WMAP result of h wo = 0.71 ± 0.04. 
As the cluster asphericity and orientation in the sky are 
random and the average cluster peculiar velocity is zero, 
these factors should contribute to the systematic uncer- 
tainty with the Hubble constant being statistically unal- 
tered. The underestimation of Hubble constant may be 
explained by the generic presence of core magnetic field 
Bo ~ 10 — 40/iG in this sample of galaxy clusters. 

Another important cosmological effect of the ubiqui- 
tous enhancement of magnetic SZE due to the prevalence 
of > l^iG magnetic fields in galaxy clusters would be ob- 
servable in the CMB angular spectrum especially at high 



I > 3000 - 4000. This contribution to CMB fluctuations 
may be estimated and tested by CMB experiments such 
as ACT, Planck, SZA, etc. Details of these two cosmo- 
logical effects will be pursued in forthcoming papers. 

For X-ray (Arnaud et al. 2001) and SZE (De Petris et 
al. 2002) spectral observations of Coma cluster (Abell 
1656), our magnetic SZE analysis is consistent with the 
currently inferred B < 10^G (Carilli & Taylor 2002). 
Likewise, magnetic SZE can be utilized in other galaxy 
clusters with high-resolution and high-sensitivity X-ray 
and SZE spectral observations to estimate the lower 
limit of Bo as well as SZE spatial features. While it 
is necessary to estimate all possible corrections to the 
classic SZE in order to isolate the magnetic contribution, 
this may be a unique procedure to probe intracluster 
magnetic fields at high redshifts z, at least statistically. 
Finally, anisotropic distributions of nonthermal electrons 
should lead to distinct magnetic SZE in radio lobes of 
extragalactic jets. 
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